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ABSTRACT 

We study recent Intermediate Luminosity Optical Transients (ILOTs) and major eruptions of Luminous Blue 
Variables (LBVs), and strengthen claims for a similar mechanism powering both. This process is a short duration 
release of gravitational energy in a binary system. In some ILOTs a merger occurs and one of the stars does not 
survive the transient event, e.g., V838 Mon and V1309 Sco. In some transient events a rapid and short mass transfer 
process takes place and the two stars survive the transient event, e.g., the Great Eruption of r\ Carinae. We study new 
ILOTs and reanalyze known ones in light of new observations and models. We reach our conclusion by analyzing these 
ILOTs using the Energy- Time Diagram (ETD) where we plot the total energy of the eruption against its eruption 
timescale. ILOTs and major LBV eruptions occupy the Optical Transient Stripe (OTS) in the ETD. The upper 
boundary of the stripe is explained by our proposed model where a main sequence (or a slightly off-main sequence) 
star accretes at a very high rate (< 1 Mq yr _1 ) from a companion. We identify one LBV, NGC 3432 OT, with two 
eruptions; one with weak total energy and the other with large total energy. It bridges the regions of the ILOTs and 
LBVs in the OTS. We further study the ILOT M85 OT2006 and show that it cannot be a nova, even not an extreme 
one. We build a model where ILOTs can become optically thin in a timescale of few years and the inflated envelope 
collapses into an accretion disk around the star. Such an ILOT will evolve blue-ward after few years. 

Subject headings: stars: winds, outflows — stars: mass loss — stars: variables: other — (stars:) binaries: general 
— (stars:) supernovae: general 



1. INTRODUCTION 

Erupting objects with luminosities between those of novae 
and supernovae (SNe) are being discovered at an accelerated 
rate (e.g., Kulkarni et al. 2007a, b; Berger et al. 2009; Smith 
et al. 2009, 2011; Kasliwal et al. 2010a; Mason et al. 2010; 
Pastorello et al. 2010). There is not yet an accepted term 
for these type of objects. We will refer to them as Intermedi- 
ate Luminosity Optical Transients (ILOTs) . Another popular 
name for these objects is luminous red novae (Kulkarni et al. 
2007a,b), although using this nomenclature may be confusing, 
as these transients are most probably not any kind of novae. 

While the nature of these eruptions is poorly understood, 
there are several different proposed explanations in the liter- 
ature (see Kashi, Frankowski & Soker 2010, hereafter KFS10, 
for further discussion). KFS10 noticed that when rescaling 
the time axis for the V-band light curves of ILOTs and ma- 
jor Luminous Blue Variable (LBV) eruptions, which at first 
seemed unconnected, the shape becomes similar for a decline 
in more than 3 magnitudes. KFS10 suggested that these 
transients may have a similar powering mechanism - mass 
accretion onto a main sequence (MS) star in a binary sys- 
tem. KFS10 also showed that the light curves are different 
than the light curves of novae and SNe (type la and II) . How- 
ever, the unique shape of the light curve of ILOTs is not yet 
understood. 
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By 'major LBV eruptions' we refer to eruptions in which 
the luminosity rapidly (< 1 month) increases by a few mag- 
nitudes, as opposed to less dramatic eruptions, such as S Dor 
phases, weak eruptions in which the luminosity is changed by 
< 0.5 magnitude in the V-band, or slow raises in magnitude. 
All these processes are most probably internally related to 
the LBVs and are not to accretion processes which we discuss 
here. 

Our proposed model of an eccentric binary interaction 
for ILOTs is directly based on the mergeburst model for 
V838 Mon (Soker & Tylenda 2003, 2006, 2007; Tylena & 
Soker 2006), in which a ~ 0.3 Mq star that merged with 
an ~ 8 Mq star created the outburst. Indirectly, the bi- 
nary model for ILOTs is supported by the similarities be- 
tween ILOTS and LBVs (KFS10), such as the massive binary 
system 77 Car. Damineli (1996) and Frew (2004; updated by 
Smith & Frew 2011) noted that the beginning of the nine- 
teenth century major eruptions of 77 Car occurred near peri- 
astron passages of the binary system. A quantitative model 
for the periastron triggering of the 77 Car nineteenth century 
major eruptions, including mass loss and mass transfer (accre- 
tion) was recently conducted by Kashi & Soker (2010). Kashi 
(2010) further argued that the eruptions of the LBV P Cygni 
in the seventeenth century were also triggered by a periastron 
passages of an invisible companion. 

The process of accretion onto a MS star was proposed in 
the past as an explanation for several different types of ob- 
jects, e.g., symbiotic stars (Kenyon & Webbink 1984). Kenyon 
& Gallagher (1985) suggested that epochs of rapid mass ex- 
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change can account for the variability of some massive stars 
although we note that r\ Car and other LBVs are consid- 
ered to be single stars in their study. We instead, suggest 
that major LBV eruptions are triggered by interaction with 
a close companion to the very massive star (Kashi & Soker 
2010; Kashi 2010). In this process, accretion of LBV mate- 
rial onto a MS (or slightly off-MS) star takes place in some 
cases in the form of an accretion disk. Bath & Pringie (1981) 
discussed eruptions from accretion episodes triggered by an 
accretion disk instability onto a MS star and a white dwarf 
(WD) in the context of dwarf novae. Though the light curve 
quite nicely matches the light curve of ILOTs when rescaling 
the time scales, the physical mechanism of ILOTs is probably 
different, as there is little chance for observing the temporary 
accretion disk through the inflated envelope. 

Smith (2011) also suggested that some (but not all) ILOTs 
are related to binary interaction and attributed the erupting 
extra energy to the motion of the companion through the pri- 
mary's envelope. It is however unclear how in such a scenario 
the energy of the eruption can be accounted for as there is 
too little orbital energy. 

In this paper we extend our model for ILOTs to account 
for recent observations. In section [2] we describe the Energy- 
Time Diagram (ETD). In section[3]we implement it for recent 
transients, demonstrating it to be a powerful tool. In section 
13] we explain why recent claims for M85 OT2006 being a nova 
are problematic, and suggest that the blue spectra observed 
in late times for M31 RV comes from an accretion disk around 
the progenitor, which is surrounded by an optically thin shell. 
We summarize in section [5] 

2. THE ENERGY-TIME DIAGRAM 

The ETD (Figure [TJ presents the total energy of the tran- 
sients, radiated plus kinetic, as a function of the duration 
of their eruptions, defined as a drop of 3 magnitudes in the 
V-band. Other energy sinks, such as the energy required to 
inflate the envelope, are not included, as they are not ob- 
served. However, in some cases the omitted energy may be 
considerable, as for the case of V838 Mon, where the energy 
required to innate the envelope is more than 10 times larger 
than the value we present in Figure [1] by the blue filled circle 
(Tylenda & Soker 2006). To account for the total available 
energy we calculate and present the gravitational energy re- 
leased by the accreted mass (marked by a black asterisk in 
the ETD). The advantage is that the total available energy 
better reflects the physical processes behind the transient, as 
according to our model these transient events are powered by 
gravitational energy. The disadvantage is that the total en- 
ergy comes from a model for each transient event, and models 
exist for only a fraction of all transient events. For LBV erup- 
tions this energy is equal to the radiated plus kinetic energy, 
as there is no inflated envelope. We estimate the total avail- 
able energy for some ILOTs as we later explain. However, 
as stated, for most ILOTs the observations and models are 
not yet detailed enough to perform this estimate, and we can 
only present the estimated radiated plus kinetic energy. For 
a pedagogical purpose and comparison, we also present in Fig 
[T]the available energy for the Type II supernova 1987A. Here 



the available energy comes from a collapsing core. Most of 
the energy is carried by neutrinos and anti-neutrinos. Only 
~ 1% of the energy goes into kinetic energy and radiation, as 
marked. 

Many transients occupy the Optical Transients Stripe 
(OTS) in the ETD, extending over ~ 3 orders of magnitude. 
The upper-right region of the OTS is occupied by major LBV 
eruptions, while the lower-left region is occupied by ILOTs. 
Figure [T] presents the updated ETD, with recent optical tran- 
sients added and other optical transient events at their up- 
dated locations, according to the values in Table [1] 

As additional transients are added to the ETD, we find it 
necessary to update the location of the OTS. The 1838 - 1858 
Great Eruption (GE) of r\ Car is presently considered as the 
highest energy and longest duration eruption in the OTS. We 
note that the 1838 - 1858 Great Eruption of r\ Car is treated 
as one major eruption, despite the fact that several peaks are 
observed throughout the duration of the outburst (KFS10). 
Smith & Frew (2011), on the other hand, treat each of the 
high rises in luminosity during the Great Eruption of r\ Car 
as a short duration transient event by itself. 

Another object, SN 2008HA, seems to be too luminous to 
belong to the OTS following the new less luminous transients 
we added (see figure 3 of KFS10). We therefore conclude that 
it most likely does not belong to this group of transient events 
and remove it from the ETD. It is important to note that the 
purpose of the OTS is not to account for peculiar types of 
SNe, and as we focus on objects which have the properties 
of accretion powered transient events we do not include such 
objects in the ETD. 

One can use a luminosity-time diagram (Kulkarni et al. 
2007a) to draw conclusions based on the observational prop- 
erties of the objects. But the usage of the ETD, with total 
energy instead of the luminosity, better accounts for the phys- 
ical properties of the different types of objects and emphasize 
that they form a stripe (the OTS) that distinguish them from 
other objects. The OTS hints to a common mechanism for 
powering transients which we propose to be associated with a 
mass accretion onto a MS star, or a star just off the MS. The 
mass transfer takes place in a binary system, although the 
companion might not survive the transient event. The OTS, 
as presented in Figure [T] has a constant slope of 1 in the 
log(£tot/ erg) - log(£/day) plane, although with the addition 
of more objects in the future the slope might change a little. 
The lower part of the OTS overlaps with novae. This should 
not be regarded as a problem for the following reasons: (1) 
Only one object, V1309 Sco occupies this overlapping region, 
and this occurs only when the energy considered is less than 
the total available energy. When all the available energy is 
considered, it is located at much higher energies and far from 
the overlapping region. (2) The ILOTs have other common 
properties that distinguish them from novae (e.g., their light 
curves) . 

The lower part is not yet theoretically constrained, and 
should be regarded as an observational constraint (selection 
effects). From above, the domain of OTS is limited by the- 
oretical considerations. Let M a and R a be the mass and 
radius of the star accreting the mass respectively (a MS or 



Table 1: Data used to update the Energy-Time Diagram (ETD; Figure[TJ. For the other objects see table 1 in KFS10. 



Transient 



Duration Total energy Total available energy Reference 

Radiated and Kinetic inch infl. envelope 
[days] [10 48 erg] [10 48 erg] 



Comments 



V1309 Sco 



PTFlOfqs 



2 r , 



NGC 3432 OT2000 62 



(5 - 15)xl(T 



> 17 0.03 - 0.3 



0.1 - 0.3 



NGC 3432 OT2008-9 531 0.85 - 2.5 

V838 Mon (updated) 70 0.03 - 0.1 

M31 RV (updated) 70 0.094 - 0.85 

M85 OT2006 (updated) 180 1.6 - 4.6 



0.023 - 0.79 

no model available 
no model available 
no model available 
0.9 - 3 

no model available 
no model available 



Mason et al. 
(2010) 



Kasliwal et al. 
(2010a); Kasli- 
wal 2010 priv. 
comm. 

Wagner et al. 
(2004) 



Pastorello et al. 
(2010) 



Tylenda & 
Soker (2006) 
Rich et al. 
(1989); Mould 
et al. (1990) 

This work; 
Kulkarni et 
al. (2007a); 
Ofek et al. 
(2008); Rau et 
al. (2007) 



a. Assuming a distance of 
8 kpc. b. Assuming the 
same ratio between total 
and radiated energy as for 
V838 Mon. 

Temporary result based on 
R-band observations. 



Light curves was corrected 
for bolometric luminosity 
using blackbody tempera- 
ture. 

Light curves was corrected 
for bolometric luminosity 
using blackbody tempera- 
ture. 



The lower (upper) en- 
ergy estimate is for 
0.01 (0.1) M Q ejected in 
the eruption. 

We estimate the ejected 
mass to be 0.2 - O.6M . 
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Energy-Time Diagram (ETD) 
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Fig. 1. — The Energy- Time Diagram (ETD). Blue empty circles represent the total (radiated plus kinetic) energy of the transients Btot as 
a function of the duration t of their eruptions (generally defined as the time it took the transient to decrease in 3 magnitudes in the V-band; 
see exceptions in the text). Blue filled circles indicate that the transient is a mergeburst. The total energy does not include the energy 
which is deposited in lifting the envelope that does not escape from the star. Where a model is available, we used is to calculate the total 
gravitational energy released by the accreted mass: the available energy. This is marked by a black asterisk. For LBV eruptions this energy 
is equal to the radiated plus kinetic energy, as there is no inflated envelope (hence the symbols overlap). In some cases, such as in V838 Mon, 
the total released gravitational energy is considerably larger than the combined radiated and kinetic energy forms. Four new eruptions that 
were added to the diagram of KFS10 are listed in Table [T] as well as updated values for eruptions which appear in that paper. Green wavy 
band: Nova models computed using luminosity and duration from Kulkarni et al. (2007a). Nova models from Yaron et al. (2005) are marked 
with red crosses. Unlike in KFS10 where we excluded part of the nova models of Yaron et al. (2005), here we present all these models. New 
models from Shara et al. (2010a) are also plotted (red open diamonds and red filled diamonds; see text for details). Eta Car GE and LE 
are the Great and Lesser Eruptions of 7] Car, respectively. The shaded area is the Optical Transient Stripe (OTS), where the eruptions we 
try to unify under the same energy source are located. The thick blue line is the upper limit of the OTS. This upper limit has a more or 
less constant luminosity as derived in equation Q. We conclude that most of these transients are associated with mass-transfer processes 
onto a MS star (we consider a merger process to be an extreme case of mass transfer in the context of our model). The LBV eruptions 
NGC 3432 OT bridge the upper and lower parts of the OTS. 
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a slightly evolved star). Star 'b' is the one that supplies 
the mass to the accretion; it can be a MS (or a little off- 
MS) star that is completely destroyed, as in the models for 
V838 Mon (Soker & Tylenda 2003; Tylenda & Soker 2006) 
and V1309 Sco (Tylenda et al. 2011), or an evolved star in 
an unstable phase of evolution that loses a huge amount of 
mass, as in the model for the Great Eruption of r\ Car (Kashi 
& Soker 2010). 

The average total gravitational power is the average accre- 
tion rate times the potential well of the accreting star 
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In the binary model discussed here, accreted mass is likely 
to form an accretion disk or an accretion belt. The accre- 
tion time must be longer than the viscosity time scale for the 
accreted mass to lose its angular momentum. According to 
Dubus et al. (2001) the viscous timescale is 
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where v is the viscosity of the disk, H is the thickness of the 
disk, C 3 is the sound speed and v$ is the Keplerian velocity. 
We scale M a and R a in equation ((2)1 according to the param- 
eters of V838 Mon (Tylenda 2005). For these parameters the 
viscous to Keplerian times ratio is \ — t v - mc /tK — 160. 

The accreted mass is determined by the details of the bi- 
nary interaction process, and differs from object to object. 
We scale it by M acc = rj a M a - Based on the modeled sys- 
tems (V838 Mon, V 1309 Sco, r\ Car) this mass fraction is 
Va i$ 0.1 with a large variation. The value of r\ a < 0.1 can 
be understood as follows. If the MS (or slightly off-MS) star 
collides with a star and tidally disrupts it, as in the model for 
V838 Mon (Soker & Tylenda 2003; Tylenda & Soker 2006), 
the destructed star is likely to be of much lower mass than 
the accretor M acc < Mb < 0.3M a . Another possibility is that 
an evolved star loses a huge amount of mass. In that case it 
is possible that the accretor will gain only a small fraction of 
the ejected mass, as in the scenario for the Great Eruption 
of rj Carinae (Kashi & Soker 2010). Here again we expect 
Afacc < 0.1M„. 

The viscous time scale gives an upper limit on the accretion 
rate 
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The maximum gravitational power is therefore 
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where we replaced the parameters of the viscous time scale 
with the ratio of viscous to Keplerian time X- Equation Q 
sets the upper bound on the OTS in the ETD, plotted as 
blue line at the upper edge of the OTS. We note that the 
location of this line may change if the accretion efficiency r] 
is different and/or the stellar parameters of the accreting star 
are different. For most of the ILOTs the accretion efficiency 
is lower, hence they are located below this line, giving rise to 
the relatively large width of the OTS. The uncertainty in r\ a 
is large and in extreme cases may be even > 1. Therefore, 
on rare occasions we expect to find objects slightly above the 
upper line drawn in the figure. 

It seems that the ETD has the potential to serve as a good 
diagnostic for classifying and perhaps better understanding 
optical transients. The next section demonstrates how we 
find recent transients to belong, or not to belong to the OTS 
and elaborate on one object that bridges the upper LBV and 
lower ILOT part. 

3. RECENT OPTICAL TRANSIENTS 

We add new objects to the original ETD of KFS10, which 
include new features that strengthen the ETD as an analyzing 
tool. We discuss each object with the new features it brings. 

V1309 Sco: The variable star V1309 Sco was discovered 
on Sep 2.5, 2008 (Nakano 2008). The distance to V1309 Sco 
was constrained to be < 8 kpc (Mason et al. 2010), and 
later better estimated as ~ 3 kpc (Tylenda et al. 2011). Ma- 
son et al. (2010) considered it to be a candidate twin of 
V838 Mon. Below we strengthen this claim for a similarity 
with V838 Mon. Tylenda et al. (2011) later conducted a thor- 
ough study of the progenitor of V1309 Sco and convincingly 
showed it to be a mergeburst similar to that of V838 Mon, and 
by that strongly supported the early claim made by Mason et 
al. (2010) and us. 

In Figure [2] we compare the V-band light curve of 
V1309 Sco (from fig. 1 of Mason et al. 2010) with that 
of V838 Mon (from fig. 5 of Sparks et al. 2008). We take 
the visible light curve of V1309 Sco (Mason et al. 2010) as 
complementary to the V-band light curve when there are no 
observations from the latter, i.e., for the first small peak at 
the beginning of the eruption. It is interesting to note that 
both objects show more than one peak at the onset of their 
eruption. We shift and scale the time scale of V1309 Sco (by 
a scale factor of 0.6) but leave the magnitude scale unchanged 
(only shifted vertically). We find that the resemblance be- 
tween the two eruptions is remarkable, as the time-scaled 
light curve of V1309 Sco follow the light curve of V838 Mon 
for a decline of almost 4 magnitudes. 

We find that the similarity between V1309 Sco and 
V838 Mon goes much further than the similarity in the shape 
of the light curves. Where data are available, the I-band is 
~ 2 magnitudes above the V-band in both objects. More- 
over, there are no high excitation lines in either object, and 
for both objects the expansion velocities are much less than 
those of novae. 

Using archive data, Tylenda et al. (2011) discovered that 
the progenitor possessed a periodic behavior with a period 
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Fig. 2. — Comparing the V-band light curves of V838 Mon 
(Sparks et al. 2008; solid-green line) and V1309 Sco (Mason et 
al. 2010). We take the visible light curve of V1309 Sco (Mason et 
al. 2010) as complementary to the V-band light curve when there 
are no observations from the latter, i.e., for the first small peak 
at the beginning of the eruption. We shift and rescale the time of 
V1309 Sco but leave the magnitude scale unchanged (only shifted 
vertically). Dashed-blue line: matching the last peak of V838 Mon 
with the second peak of V1309 Sco and then rescaling the time 
to match the light curves. One unit of time corresponds to 1 day 
for V838 Mon, and 0.6 days for V1309 Sco. The light curve of 
V838 Mon has a curved slope in the ~ 30 days following its second 
peak, before it becomes steeper. Dot-dashed red line: matching the 
initial decline of the light curves (the 'knee'). Note that the time 
scale is the same (0.6 days) for the two plots of V1309 Sco, but the 
peak magnitude and the time shift are different. Also note that the 
eruption of V838 Mon begin approximately at t = — 36 days, and 
its beginning is not shown here. Namely, the entire duration of the 
eruption of V838 Mon is longer than the eruption of V1309 Sco. 
The resemblance between the two eruptions is much better by the 
dot-dashed red line, which follow the light curve of V838 Mon for 
a decline of almost 4 magnitudes. 



of ~ 1.4 days. For the six pre-outburst years analyzed they 
found the period to decrease. They found the progenitor to 
be consistent with a binary system with a total mass of ~ 1 
3 M0 , and confirmed that the eruption was a mergeburst, sim- 
ilar to, but less energetic than V838 Mon (Soker & Tylenda 
2006, 2007). Integrating over time the bolometric luminosity 
from Tylenda et al. (2011) we find that the bolometric radi- 
ated energy is £: rad (V1309 Sco) ~ 4x 10 44 erg. This is ~ 0.017 
times the radiated energy in V838 Mon (Tylenda 2005). Since 
the bolometric radiated energy ratio between V1309 Sco and 
V838 Mon is ~ 0.02 and both are mergebursts, we take the 
total energy of V1309 Sco to also be ~ 0.02 of the total energy 
of V838 Mon. This gives £ to t(V1309 Sco) ~ 3 - 10 x 10 44 erg 
(but see the total available gravitational energy estimate be- 
low). Following this assumption we position V1309 Sco on 
the ETD (Figure [TJ. 

We note that in V838 Mon most of the energy goes into in- 
flating the huge envelope (Tylenda & Soker 2006). We there- 
fore calculate for both transients the total available energy, 
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where M a is the mass of the (more massive and dense) pri- 
mary accretor, AM ac c is the accreted mass of the compan- 
ion, and R a is the radius of accretor. We find that for 
V838 Mon the total available energy is £ to t,av a (V838 Mon) 
~ 9 - 30 x 10 47 erg, taking (M , AM acc , Raj = (8,0.3,5) in 
solar units for the lower estimate and a factor of 3.3 for the 
higher estimate. For V1309 Sco we get £tot,ava(V1309 Sco) ~ 
2.3 - 79 x 10 46 erg, taking (M„, AM acCl Ra) = (0.93, 0.07, 5.4) 
in solar units for the lower estimate, and (1.5,1.5,5.4) for the 
upper estimate, based on the results of Tylenda et al. (2011). 
As noted in section [2] the total available energy is a better 
measure of the physical process behind the accretion powered 
ILOTs. 

PTFlOfqs: The optical transient PTFlOfqs was discov- 
ered in a spiral arm of M99 by the Palomar Transient Factory 
(PTF) in Apr 16, 2010, and showed eruption characteristics 
similar to those of M85 OT2006, SN 2008S, and NGC 300OT 
(Kasliwal et al. 2010a). Taking the R-band magnitude (M. 
Kasliwal, private communication) and the spectra given in 
figure 10 of Kasliwal et al. (2010a) we find that the R-band 
luminosity should be corrected by a factor of ~ 1.1 to obtain 
the bolometric luminosity. We integrate the bolometric lumi- 
nosity over the first ~ 30 days after detection to obtain the 
total radiated energy E Iad (PTFlOfqs) ~ 3 x 10 46 erg. The 
total energy is estimated to be up to 10 times larger, about 
the maximum ratio observed for ILOTs. The peak luminosity 
occurred ~ 13 days after detection and therefore we have data 
for only ~ 17 days for the decline, during which the R-band 
decreased by only ~ 0.4 magnitudes. We note that the de- 
cline in the V-band was steeper than in the R-band. We find 
PTFlOfqs to be located on the lower edge of the OTS How- 
ever, this placement is temporary and once more observations 
are available we will be able to deduce the total energy and 
decline timescale. 

NGC 3432 OT: A series of LBV eruptions were observed 
from an LBV in the spiral galaxy NGC 3432 (Pastorello et 
al. 2010). The first in the series started in May 3, 2000 
(NGC 3432 OT2000; also referred to as SN 2000ch) and lasted 
62 days (Wagner et al. 2004). From Oct 7, 2008, a series 
of three eruptions (NGC 3432 OT2008-9) has been observed 
from the same star (Pastorello et al. 2010), lasting a total of 
~ 531 days. We refer to the above durations as the decline 
times in the ETD as they are the characteristic times of the 
transients (a similar treatment was done for the other LBVs 
in KFS10). We consider the 2000 eruption as a separate one 
from those of 2008-9 because there was a long quiescence pe- 
riod between them which had on average ~ 1.5 - 2 weaker 
V-band magnitude than at the times of the eruptions. All 
four peaks were followed by a temporary decrease in lumi- 
nosity, possibly caused by a dust shell which later dissipated, 
and allowed the remainder of the decline in the eruption lu- 
minosity to be observed (Pastorello et al. 2010). Observ- 



ing the spectra, Wagner et al. (2004) derived that the ef- 
fective temperature of the star during NGC 3432 OT2000 
was ~ 7800 K, and for NGC 3432 OT2008-9 Pastorello et 
al. (2010) obtained ~ 9000 K. More spectra given by Smith 
et al. (2011) approximately 4 years after the first outburst, 
shows Balmer lines with P Cyg profiles, blueshifted in their 
centeroids. The FWHM was ~ 1500 km s _1 with wide wings 
up to —2500 km s _1 and 4300 km s _1 in the blue and red, 
respectively (Smith et al. 2011). This is an indication of ab- 
sorbing material which was probably ejected during the LBV's 
major eruption. 

We corrected the light curves of NGC 3432 OT for 
bolomctric luminosities using blackbody temperatures for 
a star with the temperatures given above, and integrated 
the light curves. We did this procedure for the R, V, and 
B bands and derived the approximatly radiated energies of 
£ rad (NGC 3432 OT2000) ~ 10 47 erg, and £ rad (NGC 3432 OT2008- 
9) ~ 8.5 x 10 47 erg. As these transients are major LBV erup- 
tions we estimate the total energies to be ~ 3 times larger 
than the radiated energies, same as for the rj Car Great erup- 
tion (KFS10). 

NGC 3423 OT overcomes a drawback of our early analysis 
in KFS10. NGC 3423 OT2008-9 is located in the central part 
of the OTS, namely in the 2.3 < log(t/ days) < 3.3 interval 
between the upper part associated with LBVs and the lower 
part associated with ILOTs. Also, NGC 3423 OT2000 shows 
low energy and short timescale which position it at the heart 
of the ILOTs part of the OTS, though it is actually a major 
LBV eruption. NGC 3423 OT therefore bridges the two parts 
of the OTS, namely the ILOTs and the major LBV eruptions. 
This a strong support to the suggestion of KFS10 that LBV 
major eruptions and ILOTs at the lower part of the OTS have 
a similar powering mechanism. 

SN 2002bu: Smith et al. (2011) presented new observa- 
tions as well as archival study of various transients. Though 
most of them do not fall in the category of ILOTs there is one 
relevant here. 

SN 2002bu erupted in NGC 4242 (Puckett & Gauthier 
2002) and was initially thought to be a type II-P supernova. 
The V magnitude declined by 2.67 magnitude in ~ 84 days 
(Smith et al. 2011). The spectra given by Smith et al. (2011) 
shows that it reddened during the 81 days following first dis- 
covery, supposedly by a dust shell forming around it when 
the ejected material has cooled down to grain forming tem- 
peratures. This is in contrast to typical SN spectra evolution, 
leading Smith et al. (2011) to suggest that it is not a SN, 
but rather a "SN impostor". The ejection velocity, obtained 
from the Ha profile was ~ 900 km s _1 , but the ejected mass 
is hard to estimate. 

Like for V1309 Sco, we compare the V-band light curve 
of SN 2002bu from table 3 of Smith et al. (2011) to that of 
V838 Mon, and find it to follow its decline quite well when 
rescaling its time axis by a factor of 2. 

The distance estimates to NGC 4242 are not very accurate 
and we adopt the estimate of ~ 8 Mpc (Thompson et al. 
2009). accordingly, the peak luminosity translates to ~ 8.8 x 
10 40 erg s _1 . We integrated the luminosity from the time 
of discovery (~ 5 days before peak luminosity) to the last 



available observation (~ 84 days after the peak luminosity). 
Using bolometric correction to find the radiated energy, as we 
did for the ILOT in NGC 3432 has proved difficult as there 
was no conclusive blackbody temperature to associate with 
the emission. However using approximate values obtained 
from the shape of the spectra we find the radiated energy to 
be £ rad (SN 2002bu) ~ 3 x 10 48 erg. With no mass estimate 
to calculate the kinetic energy, we can only speculate that the 
total energy is in the order of ~ i0 49-50 er g ; larger than that 
of other ILOTs with the same time scale. 

As SN 2002bu is far out of the OTS, we conclude that it 
is not an ILOT, and the most likely explanation is that it 
is a peculiar class of SN. This demonstrates that the ETD 
can also be used to identify transients with different physi- 
cal properties than the ones considered as accretion powered 
ILOTs. 

4. REEXAMINING KNOWN TRANSIENTS 

4.1. Why M85 OT2006 Cannot be a Nova 

The transient M85 OT2006 was discovered in the lenticular 
galaxy M85 in Jan 7, 2006 (Kulkarni et al. 2007a). It had 
a peak luminosity of ~ 2 x 10 40 erg s _1 and total radiated 
energy of £ rad (M85 OT2006) ~ 6 x 10 46 erg over a duration 
of ~ 180 days (Kulkarni et al. 2007a). The progenitor's mass 
was estimated to be < 7 M (Ofek et al. 2008) and it was 
suggested that the origin of M85 OT2006 is a stellar merger 
(Kulkarni et al. 2007a). 

Rau et al. (2007) estimated that the effective temperature 
and the stellar radius at the time of peak luminosity were 
T e ff ~ 4600 K and R ~ 3600 Rq respectively. At a later time 
in the eruption the star cooled down and expanded to have 
Tea ~ 950 K and R ~ 20 000 R© (see their table 2). 

The work of Shara et al. (2010a) extends the previous 
nova models of Yaron et al. (2005). The data of six of 
their models is given in their tables 2-4. In the new ex- 
treme nova models of Shara et al. (2010a) the mass of the 
white dwarf is 0.4 - 0.65 Mq and the accretion rate is as low 
as 10~ 12 - 10~ 10 Mq yr" 1 . The ejected mass in the novae 
ranges between 5.3 x 10~ 4 - 2.2 x 10~ 3 M at velocities of 
150 - 480 km s" 1 . 

Shara et al. (2010a) claim that the eruption of M85 OT2006 
was a nova, and suggest that it can be explained in the frame 
of the new extreme nova models. Their conclusion is mainly 
based on the following. (1) Their new results showing that 
novae can reach comparable luminosities of ILOTs such as 
M85 OT2006 and M31 RV (few xlO 7 L ). (2) The models 
produce red eruptions, as the spectra of ILOTs. 

We hereby show that the mass ejected in the eruption of 
M85 OT2006 is much larger than the nova models can pro- 
duce. Let us consider the M85 OT2006 transient according 
to data given by Rau et al. (2007), discussed above. The 
column density required above the photosphere is given by 

AA~pAr = \- =0.67 g cm" 2 (- \ T ) (6) 

3 k \ 1 cm 2 g / 

where p is the average density above the photosphere, Ar is 




Fig. 3. — The photospheric mass and its opacity expected 
from M85 OT2006 at (or right after) peak luminosity, when ef- 
fective temperature is T ef j ~ 4600 K and the stellar radius is 
R ~ 3600 Rq, plotted as a function of the relative thickness of 
the photosphere /3 = Ar/R. Different cases are for different values 
of Hydrogen abundance X and metallicity Z (see legend) . 

the thickness of the shell above the photosphere, and K is its 
average opacity. 

We parameterize the thickness of the shell with Ar = fiR 
with p ~ 0.1. The final mass we obtain for the shell above the 
photosphere changes by a factor < 2 for 0.01 ^ ft ^ 0.1. The 
reason is because the value of ft determines the density that 
weakly influences the opacity. We use opacities from Ferguson 
et al. (2005) (slightly extrapolated), using compositions from 
Asplund et al. (2004) with hydrogen abundance X = 0.7 and 
metallicity Z = 0.1; other composition from Lodders (2003) 
gives very close opacity values. The total mass above the 
photosphere is 

M p h = 4irR 2 pAr 

^xio-r^-Vr : r v x M S . (7) 

Vl0 4 R e J ^lcmag-V ° 

Using the temperature at peak luminosity we find that the 
opacity at the peak luminosity is k ~ 1.3 x 10 -3 . Substitut- 
ing the opacity and the radius at peak luminosity in equation 
(0 we find that the mass above the photosphere at the time 
of peak luminosity is M p h ~ 0.2 Mq. Our results are sum- 
marized in Table [2] The real amount of mass is even larger, 
as some mass is well above the photosphere and we expect 
a large amount of mass to be below the photosphere as well. 
The largest value of ejected mass in the nova models of Shara 
et al. (2010a) is 2.2 x 10~ 3 M or two orders of magnitude 
below the expected ejected mass we calculate. We therefore 
conclude that M85 OT2006 is not a nova. 

Previously in KFS10, we used the assumption of Ofek et al. 
(2008) that the total ejected mass of M85 OT2006 is 0.1 M , 
and consequently obtained a total energy of ~ 1.4 x 10 47 erg. 
Our new estimate of ejected mass is much higher and we up- 
date our estimate of the total energy in M85 OT2006. Our 



new estimate, taking ejected mass velocity of ~ 870 km s 
(Rau et al. 2007) is 1.6 x 10 48 - 4.6 x 10 48 erg, corresponding 
to 0.2 Mq and 0.6 Mq, respectively. These values are sum- 
marized in Table [1] and the updated location of M85 OT2006 
is shown in the ETD (Figure [l} . 

We calculate the total (radiated and kinetic) energy of the 
new nova models of Shara el al. (2010a) and plot them in 
Figure [T] We include all models from Yaron et al. (2005) (red 
crosses) including those we omitted in KFS10. The data of 
Shara et al. (2010a) is missing the parameter t m i, the duration 
of the mass-loss phase, which was given in the previous list of 
Yaron et al. (2005). The six extreme nova models of Shara 
et al. (2010a) are therefore plotted twice in Figure [T] First, 
taking i3,boi instead of t m i in the calculation for the total 
energy (red open diamonds). Second, following rough average 
values from Yaron et al. (2005), that show the relation t m \ ~ 
0.1£3,boi (red filled diamonds). We find the new extreme nova 
models to be located in the lower, less populated part of the 
OTS. This implies that novae are very unlikely to account for 
the ILOTs. 

4.2. M31 RV: Seeing an Accretion Disk Through 
the Ejecta 

The outburst of M31 RV occurred in 1988 and had a 
peak luminosity of ~ 8 x 10 Lq, declining on a timescale 
of ~ 70 days (Rich et al. 1989; Mould et al. 1990). The 
progenitor is located in the central bulge of M31 (Rich et 
al. 1989), where the stellar population is old. The ejecta 
was estimated to have a mass of 0.001 - 0.1 Mq and an ex- 
pansion velocity of 100 - 500 km s" 1 (Mould et al. 1990). 
About 54 days after the peak the spectra of the ILOT was 
compatible with a blackbody temperature of ~ 2000 K, and 
79 days after the peak it dropped to ~ 1050 K (Mould et al. 
1990). Radiation from a ~ 1000 K dust shell with a radius of 
~ 8000 Rq was observed ~ 70 days after the eruption, infer- 
ring that the average expansion velocity over this period was 
~ 920 km s" 1 . 

Based on some similarities with V838 Mon, Soker & 
Tylenda (2003; Tylenda & Soker 2006) suggested that 
M31 RV eruption was caused by mergeburst. According 
to the theme of the present paper, a merging processes can 
be replaced with a rapid mass transfer episode in a binary 
system. In a recent paper, Shara et al. (2010b) analyzed HST 
archival data, and found that in 1994 the star's luminosity 
was ~ 10 3 Lq and its temperature was > 40 000 K. Shara 
et al. (2010b) also presented new HST observations of the 
surroundings of M31 RV from 2008, and claimed to detect 
an energy source with an effective temperature of ~ 8000 K 
at the location of the progenitor. According to Shara et 
al. (2010b), there is currently no nova model that can pro- 
duce the observed peak luminosity of M31 RV. However, 
they suggested that with a usage of the correct opacity (that 
presently their model lacks) low mass white dwarfs accreting 
at low rates can in principle produce low temperature and 
very luminous novae, as required in the case of M31 RV. 

Motivated by the new HST observations of Shara et al. 
(2010b) we consider transients with low ejected mass. Lo- 
cated among old stars, the mass of M31 RV is expected to be 



Table 2: Estimate of the mass above the photosphere in the eruption of M85 OT2006. Data for effective temperature and photospheric 
radius are taken from Rau et al. (2007). Opacities are from Ferguson et al. (2005) (slightly extrapolated), using the composition from 
Asplund et al. (2004) with Hydrogen abundance X = 0.7 and metallicity Z = 0.1. The estimated mass above the photosphere is calculated 
from equation JTJ • Results are given for two values of the relative thickness of the photosphere f3 = Ar/R. 



Time 


T eff [ K] 


R{ R ] 


k[ cm 2 g 1 ] 


M ph [M ] 








(3 = 0.1 p=l 


13 = 0.1 /3 = 1 


Peak 


4600 


3600 


1.3 x 10- 6 1.0 x 10-* 


0.20 0.26 


Late 


950 


20 000 


4.45 2.35 


1.8 x 10" 3 3.5 x 10" 3 



low. Possibly, the ILOT is a merger event of two low mass 
stars or a mass transfer episode to a low mass star. In such 
cases the mass transferred to the accreting star is low, and the 
ejected mass is even lower. We suggest that the ejected mass 
was small, and quickly became optically thin. Part of the 
ejected material might have not reached the escape velocity, 
and fell back toward the star. 

In addition, in our model the inflated envelope has also 
contained very little mass (< 0.01 Mq). This is different than 
the situation in V838 Mon, where the inflated envelope was 
more massive, and is long-lived (Tylenda, Kamiski & Schmidt 
2009). In the case of V838 Mon ~ 0.1 - 0.3 M has been 
accreted on a radiative envelope of a B-star (Tylenda, Soker 
& Szczerba 2005; Tylenda & Soker 2006). On the other hand, 
in the case of M31 RV, the stellar envelope is convective, and 
as we show below that only ~ 0.05 Mq has been accreted 
onto it. The convective envelope has quickly adjusted itself 
to the accreted mass. A large envelope with a very low mass 
cannot support itself (as is known for AGB stars for example; 
Soker 1992), and therefore the inflated envelope in the case of 
M31 RV has collapsed onto the low mass star. 

The gas in the inflated envelope that later collapsed, and 
the gas that did not reach the escape velocity at larger dis- 
tances, fell back toward the star. In case of a surviving com- 
panion, the companion could have continued to transfer mass 
to the mass accreting star. In any case, instead of a long-lived 
inflated envelope as in the case of V838 Mon, an accretion disk 
was formed after about a year. 

We can propose some typical numbers, but with very large 
uncertainties. The fall back mass was Mfb ~ 10 -3 M . To 
account for a luminosity of ~ 1000 L after several years 
(Shara et al. 2010b) an accretion rate of ~ 10 -5 M yr -1 is 
required. The disk can live for ~ 100 years. 

Because of the high specific angular momentum in the bi- 
nary system the fall back gas formed an accretion disk around 
the accreting star. After a relatively short time (~ 1 yr in 
the case of M31 RV) one would observe the central star and 
its accretion disk rather than a huge inflated envelope. The 
accretion disk and its boundary layer become the dominant 
illuminating source. 

We start by following the calculation of Shore & King 
(1986), who studied the physics of accretion disk boundary 
layers. For a star with mass M a and radius R a , accreting at 
a rate of M a , the luminosity of the boundary layer would be 
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Having an optically thick boundary layer with a thickness 
<5bl, it will heat to a temperature of 
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(9) 



where a is the Stefan-Boltzmann constant. 

The temperature of the boundary layer for the parameters 
used in equations |8} and Q is compatible with the hot tem- 
peratures of M31 RV as seen in the observations of Shara et 
al. (2010b) at late times. The accretion luminosity is evenly 
divided between the boundary layer and the accretion disk 
itself. The accretion disk temperature is much lower than 
the temperature of the boundary layer. Hence, it is expected 
that the averaged observed temperature will be lower than 
the temperature of the boundary layer. 

Let us now check what would be the optical depth of the 
transient a few years after its eruption. We will assume spher- 
ical symmetry. The thickness of the shell, Ar er , ejected during 
the eruption is determined from the timescale of the eruption 
At cr in which the material is ejected, and the velocity v aT it 
had during the eruption 



Ar er ^ 20 ( — - 

\ 500 kins- 1 



Ater 

70 days 



AU. 



(10) 



We take the mass of the spherically expanding shell M e j, and 
its average expansion velocity w e j. After a time period At it 
would reach a distance of 



^ 1200 1 920 km S - 



At 
6 yr 



AU, 



(11) 



where the calibration is from M31 RV and the latest HST 
observations of Shara et al. (2010b). The average density in 



the shell is therefore 



Me 



47rr|jAr e 
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20 AU 
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and the column density of the shell is 



AN c j ~ PejAr e 



5.2 x 10" 



M c 



0.01 M 



1200 AU 



g cm 



(13) 



For a shell temperature of 1000 K with Hydrogen abun- 
dance X = 0.7 and metallicity Z = 0.02 the opacity 
(from the same data by Ferguson described above) is K c j = 
0.03 cm 2 g _1 , and the optical depth is 



1.6 x 10" 4 



AN C : 



5.2 x 10" 3 g cm- 2 J V0.03 cm 2 g" 1 



(14) 



For temperatures in the range ~ 600 - ~ 3000 K the opac- 
ity is K j < 1 and therefore the result that the shell be- 
ing optically thin is not sensitive to temperature. If dust 
is formed the opacity can reach much higher values. Several 
researchers have suggested that dust opacities can reach val- 
ues of ~ 10 cm 2 g _1 and even higher (e.g. Helling et al. 2000; 
Henning & Stognienko 1996; Pollack & Mckay 1985; Pollack 
et al. 1994; Semenov et al. 2003). However, these high opac- 
ities were not obtained for densities as low as ~ 10 -17 , as in 
our case. According to the studies discussed above for such 
densities, the opacity is K e j < 1 in the temperature range 
~ 100 - ~ 1000 K, which is the relevant range for our study. 
From equations (fTI)) . (fT2"1) . (fT3)l and (fTI)) . the time period it 
would take for the dusty shell to become optically thin is 



1 / K C \M C \ 

At = — ' J J 



0.4 | 



920 km s- 



1 cm 2 g - 



(15) 



M c 



0.01 M. 



— ) 



We see that even if dust grains are formed, a few years af- 
ter the eruption r e j << 1, and the shell becomes optically 
thin. This result holds even if the ejected mass is as high as 

A'/ cj , m ax = 0.1 Mq. 

We therefore propose the following scenario for the out- 
burst of M31 RV. This scenario with different parameters is 
relevant to other ILOTs with low mass progenitors as well. 
An interaction between two old stars, a primary with a mass 
of ~ 1 Mq and a less massive companion, led to the ejection 
of material that created a shell around the primary. Alterna- 
tively, this interaction might have been a mass loss episode 
of the companion due to tidal interaction with the primary. 
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Approximately, a mass of M e j ~ 0.01 Mq was ejected and cre- 
ated an expanding shell around the stars. The energy to eject 
the shell comes from accretion of material onto the star(s). 
The kinetic energy of the ejecta is 



E elMn ~ 8.4 x 10 4 



M c 



:)( 



0.01 M 7 V920 km s" 



erg. (16) 



The radiated energy £ cr , rad ^ 10 46 erg (Mould et al. 1990). 
Therefore the total energy of the eruption is -E cr ,tot = £'er,kin+ 
£ er , rad ~ 10 47 erg. We assume that the companion which 
donated the accreted and ejected material had a highly ec- 
centric orbit e > 0.85, and that the accretion episode oc- 
curred close to periastron where the tidal force was maxi- 
mal. Therefore, the accreted material fell onto the star at 
approximately the free-fall velocity. Under this assumption 
and according to the virial theorem, the potential energy was 
(U) — — 2(£ C r,tot) = 1.9 xlO 47 erg. The accreted mass needed 
to supply this amount of energy is 



M acc 



2-Eer.tot Ra 



GM a 



0.05 



E c 



10 X 10 47 ere 



-Ra 

1 R e 



1 M(7 



Ma 



(17) 



In other words, most of the gravitational energy of the ac- 
creted material went to kinetic energy to eject the shell, and 
a small part went to radiated energy. 

We conclude that after a time period At ~ 1 yr the shell 
became optically thin and the with lack of obscuration from 
an inflated envelope the disk became observable (equation 
I15|l . Therefore, a few years after the eruption a blue energy 
source is seen at the location of the progenitor, as observed 
by Shara et al. (2010b). 

5. SUMMARY 

In a previous paper (KFS10) we grouped several Interme- 
diate Luminosity Optical Transients (ILOTs), together with 
major Luminous Blue Variable (LBV) eruptions. These ob- 
jects occupy the Optical Transients Stripe (OTS), between 
novae and supernovae in the Energy- Time Diagram (ETD; 
Figure [JJ, where the transients' total energy (kinetic and ra- 
diated; or total available energy when a model exist) is plotted 
against their decline timescale. 

Recent observations of new objects (section fallow a bet- 
ter placement the OTS on the Energy-Time plane. The LBV 
eruptions of NGC 3432 OT2 (Pastorello et al. 2010), bridge 
the gap between the upper and lower parts of the OTS (see 
Figure [TJ. NGC 3432 OT serves as strong support for group- 
ing LBV major eruptions together with ILOTs. This allows 
us to overcome a drawback of the original proposal, hence 
strengthening the usage of the ETD. The upper bound of the 
OTS comes from accretion considerations and has a constant 
luminosity of ~ 10 42 erg s _1 , for a maximal accretion rate of 
~ 10% of the accretor mass (equation [4]). The lower part of 
the OTS is only observationally constraint, and has no the- 
oretical limit. It is however evident that the luminosity of 



ILOTs is larger than novae despite a small overlap (which we 
explain in section [2]|. 

An important contribution to our physical understanding 
of ILOTs is the new ILOT V1309 Sco (Nakano 2008) and its 
similarity to the eruption of V838 Mon, as noted by Mason et 
al. (2010) and later confirmed by Tylenda et al. (2011). The 
strong similarity in light curves (Figure [2] after time rescal- 
ing) hints to a basic common physical mechanism. We fully 
accept the mergeburst model of V838 Mon and V1309 Sco, 
and we are not proposing an alternative. On the contrary, we 
extend the mergeburst model by considering a massive and 
short duration mass transfer episode that do not destroy the 
mass losing star. A merger process is considered to be an 
extreme case of mass transfer events. 

We refute claims by Shara et al. (2010a) that the ILOT 
M85 OT2006 was a nova in light of recent extreme nova mod- 
els. We show that the energy of M85 OT2006 is much too high 
for any nova (even extreme nova), and place it in the ILOTs 
part of the ETD (section |4T]). Another ILOT, M31 RV, was 
observed in an old stellar population. This implies that the 
processes which lead to its eruption involved low mass stars. 
In our accretion model, a relatively small amount of mass 
was ejected from the system, and a small amount of mass was 
in the inflated envelope. The small amount of ejected mass 
ensured that the shell became transparent ~ 1 yr after erup- 
tion. In addition, with little mass in the inflated envelope, 
the envelope collapsed within few years. The fall-back mate- 
rial in our binary model has high specific angular momentum, 
such that an accretion disk is formed around the star. The 
accretion disk and its boundary layer dominate the radiation. 
The temperature is relatively high, ~ 5000 - 50 000 K. This 
explains the blue-ward evolution of M31 RV years after erup- 
tions (section |4. 2 [) . as observed by Shara et al. (2010b). 

More transient objects will very likely be observed in com- 
ing years. The OTS is a strong tool in analyzing these tran- 
sients, as well as light curve matching as we demonstrate here. 
Other research groups have also begun to adopt energy consid- 
erations. For example, the .la supernova SN2010X (Kasliwal 
et al. 2010b) was analyzed in terms of its total energy, rather 
than its peak magnitude, as we suggested in KFS10. We en- 
courage other groups to do the same, and to check whether 
the transient events they observe are located on the OTS. 
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